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FIXED POINT FOUR-FERMI THEORIES 


Simon HANDS 

Department of Physics, University of Wales Swansea, 

Singleton Park, Swansea SA2 8PP, UK 

I review dynamical chiral symmetry breaking in four-fermi models, including results of 
Monte Carlo simulations with dynamical fermions. For 2 < d < 4, where the phase 
transition defines an ultraviolet fixed point of the renormalisation group, the continuum 
theory may either be describable using the large-Aj^ expansion, as in the case of the 
Gross-Neveu model, or be intrinsically non-perturbative, as in the case of the Thirring 
model. For d = 4, the models are trivial and are described by a mean field equation of 
state with logarithmic corrections to scaling, which may nonetheless define new univer¬ 
sality classes distinct from those of ferromagnetism. 


1 Introduction 

In this talk I will review recent progress in our understanding of four-fermi theories, 
by which I mean quantum field theories described by the Lagrangian density 

r denotes a 4 x 4 matrix constructed from Dirac 7 -matrices, and the sign of the 
interaction is chosen so as to be attractive between fermion and anti-fermion. The 
index i runs over Nf distinct fermion flavors. I will discuss these models in spacetime 
dimension 2 < d < 4, paying most attention to the cases d = 3 and d = 4. 

What are the motivations for studying such models? Their generic behaviour 
is that dynamical breaking of a chiral symmetry occurs at some strong value of 
the coupling ( 7 ^, where non-perturbative methods must be applied. In this talk I 
shall discuss and compare three such methods: the large-V/ expansion, use of the 
Schwinger-Dyson Muation, and lattice Montje Carlo simulation. Other methods, 
such as variationallH and derivative expansioncl approaches, have also been applied. 
The picture which is emerging is that the symmetry breaking transitions in these 
fermionic models define new universality classes, resembling qualitatively, but not 
quantitatively, those which apply to the study of ferromagnetic phase transitions. 

Phenomenologically, four-fermi theori^ were of course originally introduced 
by Fermi to describe the weak interaction 0, and were next applied to dynamical 
chiral symmetry breaking in the strong interaction by Nambu and Jona-Lasiniou. 
More recently they have appeared in discussions of dynamical mass generation in 
the Standard Model, in such scenarios as walking technicolom, and the top-mode 
standard model in which the Higgs scalar is a tt bound state 0. In three spacetime 
dimensions, there may be applications to high-T^ superconductivity!], for instance 
in describing non-Fermi liquid behaviour in the normal phase □. 
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2 The Gross-Neveu Model for d = 3 

I will begin by discussing the simplest model, the Gross-Neveu (GN) model i, in 
which most of the important theoretical issues are already present. The Lagrangian 
is ^ 

Cgn = '0i(^ + rn)tlji - ■ (2) 

For bare fermion mass m = 0, there is a discrete chiral symmetry 


tjj 75 ?/) ; ij} -?/) 75 , 


( 3 ) 


which is spontaneously broken whenever a non-vanishing condensate (t/'V’) is gen¬ 
erated. To proceed, we introduce a bosonic scalar auxiliary field a and rewrite 

Cgn +m + a)'ijji + (4) 

The original Lagrangian (H) can be recovered by Gaussian functional integration 
over a. Chiral symmetry breaking for m ^ 0 is now signalled by a non-vanishing 
vacuum expectation value S = (cr) for the scalar field: it then follows from (^) that 
the fermion gets a dynamically generated mass M ~ E. 

We can calculate E using an expansion in inverse powers of fV/, the number of 
flavors. This expansion associates a factor Nf with each closed fermion loop, and 
in effect l/^/Nf with each fermion-scalar interaction vertex. To leading order, in 
the chiral limit m 0, only the tadpole diagram shown in Fig. la contributes to 
E, leading to the self-consistent Gap Equation: 


E f 1 


9 


( 5 ) 


or, with a simple UV cutoff A on momentum (note 2 < d < 4) 


r (4^)i(d-2) 


Ad-2 

nf) 


_Ed-2r(2_ I) 


( 6 ) 


Equation (||), relating a bare coupling constant to both a UV scale A and 
a physical scale E, can be interpreted as a renormalisation condition. It turns out 
to be the physical solution of (||) for couplings larger than the critical coupling 
given by 

^ ^ 8 A^-^ 

5? (47r)f(d-2)r(|)’ 

at which point chiral symmetry is spontaneously broken (see Fig. 4). Only for 
g^ —> g^, can the ratio A/E be made to diverge, implying a continuum limit. Note 
that it is also possible to approach the continuum limit from the symmetric phaseE3. 

Remarkably, for 2 < d < 4, to the same leading order in 1/Nf all dependence on 
A is absorbed in defining the value of 5 /. The only other Green function involving 
a fermion loop is the scalar two-point function shown in Fig. lb, but in the vicinity 
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(a) 


(b) 


Figure 1: Leading order diagrams in the GN model 


of the fixed point, once (|^) is taken into account the scalar propagator is finite and 
can be expressed in closed formEdl. Eg, for d = 3, 


D.{k) = 


1 2TT^yk^ 

+ AY?) tan“^ (’^^) 


( 8 ) 


In the large-A^/ limit is essentially a chain of fermion bubbles. It is worth 
examining its behaviour in two limits. In the infra-red, 


limZIcr(fc) oc 
fc-»0 


1 

A:2 + 4E2 ’ 


(9) 


and hence the a resembles a fundamental boson with mass 2Y. This shows the 
scalar to be a weakly bound fermion - anti-fermion composite. In the ultra-violet, 
on the other hand, for general d we have 

oc (10) 


Thus the UV asymptotic behaviour is harder than that of a fundamental scalar, but 
still softer than the l//c° corresponding to a non-propagating auxiliary field. The 
strong interaction between fermion and anti-fermion is responsible for this modi¬ 
fication, which in turn makes diagrams corresponding to higher order corrections, 
such as those of Fig. 2, less divergent than expected by naive power counting. 

Thn result, known for a long time 113 , but discussed recently with renewed in¬ 
terest E3l, is that the 1/Nf expansion about the fixed point = g1 is exactly 
renormalisable for 2 < d < 4. 

An important aspect of the model’s renormalisability is that it depends on a 
precise cancellation between logarithmic divergences from different graphsllil. This 
can be understood in terms of the four-fermi scattering amplitude, shown schemat¬ 
ically in Fig. 3: there are three different types of logarithmic divergence, each 
represented by a blob, but in the massless limit only two tuneable parameters in 
the renormalised Lagrangian, implying a non-trivial consistency relation. 

In turn this implies that in the deep UV limit the amplitude assumes a universal 
form 


Jim Mff^ff 

k —»-oo 


Ad 

Nfk^-"^ ’ 


( 11 ) 
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Figure 2: 0(1/Nf) radiative corrections in the GN model 





Figure 3: Renormalised four-fermi scattering amplitude 


with the numerical constant Ad = /4r(2 — — 1), and crucially no 

l/Nf corrections. How then do higher order corrections manifest themselves? Fig. 4, 
showing a comparison between Monte Carlo data for S from a 12^ lattice with the 
predictions of the leading order gap equation (||), gives a clue. The corrections to 
E are clearly 0{1/Nf), suggesting that it is the behaviour of the bulk “thermody¬ 
namic” observables near the phase transition, in other words the IR physics, which 
most conveniently expresses the higher order information. If we define a reduced 
coupling 
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( 12 ) 


then the scaling of the order parametet,^near the fixed point may be described by 
a set of critical exponents. Explicitly tilEJ: 


Soct^ with f3=^ + o(^-^y, 


E|t=o oc with 6 = {d — 1) 


1-b 


Nf 


X = 


as 

dm 


oc \t\ with 7 = 1 - 1 - 


m—0 


d-lCd, 
d-2 Nf’ 


M oc 1^1“^ with n = 


1 


d-2 


1-b 


d-lCd 


d Nf\ 


(13) 

(14) 

(15) 

(16) 
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Figure 4: Plot of E vs. for Nf = 24 (squares), 12 (circles), and 6 (triangles). The line is the 
leading order gap equation prediction. 


n^(fc)|t=o ocwith 77 = 4-d-(17) 

The numerical constant Cd = — 1). Several comments can be 

made: 

• The exponents obey certain consistency checks known as scaling and hyper¬ 
scaling relations, eg. 

7 = 1 /( 2 - 77 ) ; /3 = f 7 /(d -2 + 77 ). (18) 

In statistical mechanics these are derivedlli on the assumption that there is a 
single length scale characterising the important physics, namely the correla¬ 
tion-length With ‘C ~ this is precisely the statement of renormalisabil- 
ityld. 

• In the limit d ^ 4 the exponents assume their mean field (Landau-Ginzburg) 
values, namely 

/3=i ; (5 = 3 ; 7 = 1 ; 7 /=^ ; 77 = 0 . (19) 

• The exponent 77 is relcded to the anomalous scaling dimension of the 
composite operator 

f] = d-2y^^. ( 20 ) 
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Table 1: Critical exponents for the d = 3 Nf = 

2 Gross-Neveu model 

rn 


large-A/EfI 

£-expansionEi3 

Monte Carlolli^ 

V 

ijv 

0.903 

1.2559 

0.9480 

1.237 

1.00(4) 

1.246(8) 


This implies for the scaling dimension the following: 

= d-1-(d-2 + 0(l/Af/)) 
= 1 + 0(1/7V/). 


( 21 ) 


Hence the interaction term has scaling dimension ~ 2, and is super- 

renormalisaiJe, or in the renormalisation group (RG) sense relevant at the 
fixed point E3.|J^ote that the sign of the 0{1/Nf ) correction is not crucial to 
the argumentEiJ. 

With the addition of and cr'* terms, the model becomes a Higgs- 

Yukawa (HY) theory, which is super-renormalisable for d < 4. The UV fixed 
point of the GN model can thus be viewed as an IR fixed-point of the HY 
model, at which the two extra operators become irrelevant til. This recalls the 
Wilson-Fisher fixed point in ferromagnetic models, and suggests an expansktn 
in £ = 4 — d as an alternative approach to the study of the universality clas^ld. 
The equivalence of the two models in d = 3 has been tested numericallyEj. 


The final issue we should address for the GN model is the applicability of the 
1/7V/ expansion; in other words, for how small a value of Nf may we trust its 
accuracy? Calculations of critical exponents are now available to 0(1/Y'|), and are 
compared to predictions from the e-expansion of the HY model to O(e^), and the 
results of a d = 3 Monte Carlo simulation for = 2 (the smallest value of Nf that is 
simulable with a local lattice action using the staggered fermion formulation, which 
retains a chiral symmetrylll^) in Tab. |^. The agreement is satisfactory, though it 
should be noted that the 1/iV/ expansion for v, and the e-expansion for ^ jv appear 
slowly if at all convergent, anc^he values quoted in the table have been obtained 
using resummation techniques^. 

A similar comparison between the large-A^/ expansion to 0(1/Af|) and Monte 
Carlo results has also been made for the GN model with continuous U(l) chiral 
symmetry, for Nf = with the results shown in Tab. |^. In this case no resummation 
is attempted on the 1/A^/ series: instead the quoted error reflects the size of the 
last term. 

We conclude that the physical picture of dynamical symmetry breaking in the 
GN model given by the large-expansion is borne out by explicit simulations 
in d = 3, and is qualitatively correct. Quantitative accuracy is limited by the 

Th^ formulation has a parity-invariant mass term, consistent with our use of four-component 
spinors Ell. 
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Table 2: Critical exponents for the d = 3 Nf = 4 Gross-Neveu model with U(l) chiral symmetry 



large- A^/E j 

Monte Carlo0 

u 

1 . 0 ( 1 ) 

1 . 02 ( 8 ) 

7 / 1 / 

1.055(13) 

1.19(13) 

Plu 

0.973(6) 

0.89(10) 


apparently rather slow convergence of certain series, such as that for the exponent 

V. 

3 The Thirring Model for d = 3 

The next model to consider, which will turn out to be more interesting, ja the 
Thirring model, originally studied as a soluble model for d = 2 and Nf = 1 E3, but 
here generalised. It differs from the GN model in that the contact interaction is 
between conserved currents: 


jO-Thir = 

INf 

Nf 

= '0i(^ + *4 + (22) 

where in the second line we have introduced an auxiliary vector field A^. The 
form of the first term of ( p^ suggests that of an abelian gauge theory, except that 
the term quadratic in violates the gauge symmetry. In the chiral limit the 
Lagrangian ( p^ ) has a continuous U(l) chiral symmetry: 

Ip exp(iaj5)ip ; tp ip e'xp{iaj5) . (23) 


It is possible to develop the 1/iV/ expansion just as we did for the GN model, but this 
time in the chiral limit (ipip) remains zero to all orders, essentially because the trace 
over an odd number of ^matrices vanishes. Gonsequently there is no fixed-point 
condition analogous to (g); however, provided a regularisation is specified which 
preserves current conservation d^{'ipjfj,tp) = 0 , the auxiliary propagator D^i,{k) is 
actually finite at leading order. For d = 3 we have, explicitly 


= ( (5^. - 


5" 


k^ki^ 


I + — \ m + 

ZTT 


k^ — , / Jk'^ 

tan-i ' ^ 




V 2m 


k^ki, 
fc2 ■ 


(24) 

The longitudinal piece has no physical content, since it vanishes on being applied 
to a conserved current. As before, we can examine (M) in two opposite limits. In 
the UV, we find 

1 


lim D^v(k) oc ( 


^2 j J^d-2 ' 


(25) 
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Figure 5: Schwinger-Dyson equation for the full fermion propagator 


Once again, the modified asymptotic behaviour means that the divergence structure 
of higher order corrections is better behaved-thap-expected: the Thirring model is 
thus exactly renormalisable for 2 < d < In the IR limit, the position 

of the pole in the transverse piece of yielcls_.the mass My of the propagating 
boson in the vector channel. For d = 3 we havec3 


Mv 


2m (1 — 2 exp 




^ 0 ; 

mg'^ oo. 


(26) 


The result depends in a non-trivial way on the dimensionless combination mg'^. In 
the weak coupling limit mg^ —> 0 the boson is a weakly bound fermion - anti-fermion 
state as in the GN model; however in the strong coupling limit its mass vanishes, 
yielding a model in which a massless vector mediates interactions between conserved 
currents. This is strongly suggestive that in this limit the Thirring model is identical 
to QEDg. Indeed, in general we expect the UV limit of the Thirring model to 
coinicide with the IR liipi^f QEDg, since in these limits the 1/Nf expansions of 
the two models coincideB3E3. 

We can summarise the distinct features of the Thirring model in 2 < d < 4, 
as compared to the generic GN behaviour described in the previous section, with a 
number of equivalent statements: 


• There is no fixed-point condition. 

• The bare parameter g^ is not renormalised. 

• The continuum theory contains a free dimensionless parameter mg'^. 

• The interaction is marginal in the RG sense. 

Is this the whole story? All of these statements are made plausible in the 
context of the large-iV/ limit. We will now argue that the behaviour of the Thirring 
model is radically different for small JVf, if methods which are non-perturbative in 
Aif are used. ^ ^ I—. 

The Schwinger-Dyson (SD) approach cZl e 2I 1^ is to solve for the full inverse 
fermion propagator = A(k)(i^ +M(k)) self-consistently using integral equa¬ 
tions, shown diagramatically in Fig. 5. 

Any non-trivial solutkin M (k) ^ 0 in the chiral limit signals dynamical mass 
generation. One analysisexploits a hidden local symmetry of to fix a non¬ 
local “gauge” in which the wavefunction renormalisation A{k) = I. To proceed, the 

^Stpeitly renormalisablity beyond leading order has only been demonstrated for the massless 
theory ESi. 
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Figure 6: Proposed UV fixed points for the small-Thirring model 


SD equations must be truncated; the full vector propagator is set to its large-Nf 
form and the full vertex r^(p, q) to the bare vertex - this second step is 
known as the ladder approximation, since in effect it ignores diagrams with crossed 
vector lines. With these approximations the SD equation may be solved in closed 
form in the strong coupling limit —> oo. 

The result is that a non-trivial solution is found for small values of Nf] explicitly, 
for d = 3 


M 

— oc exp 



^ cx A^exp 



(27) 


with A a UV cutoff. The solution exists for 

1 28 

4.3. (28) 

If Nf were a continuous parameter, then ( p^ would predict a UV fixed point 
as Nf^ ^fc- The essentially singular behaviour is similar to that of quenched 
QUD , Pd; this type of transition has been given the label “conformal phas e transi¬ 
tion”!^. How are we to interpret it in this case? It has been argued Edo that for 
< oo the critical curve Nfc{g) is smooth. It is tempting to propose the phase di¬ 
agram shown in Fig. 6, implying a sequence of fixed points g‘^{Nf) for integer values 
of Nf < Nfc- The Thirring interaction (Uy^i'i/’)^ has become relevant. We-should be 
cautious, however; a different sequence of truncations predicts Nfc = ooEll. Indeed, 
the situation is reminiscent of QEDg, where chiral symmetry breaking is predicted 
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Table 3: Fits to the Equation of State 



Nf=2 

Nf=4. 

Nf = 6 

Wc 

1.92(2) 

0 . 66 ( 1 ) 

no fit 

6 

2.75(9) 

3.43(9) 

found 


for Nf < Nfc, the predicted valufi-of Nfc being highly sensitive to the assumptions 
used in solving the SD equation E3l3. Once again, however, the same SD equation 
appears to describe both models, but in opposite momentum regimes, suggesting 
that their equivalence may hold beyond the large-approximation. 

The UV behaviour of the Thirring model, if non-trivial, is genuinely non- 
perturbative, since there is no small dimensionless parameter in play at the pu¬ 
tative fixed points. .In, this sense the Thirring model resembles the strong-coupling 
behaviour of QED 4 E®lEZl, and hence is a suitable case for lattice Monte Carlo treat¬ 
ment. Indeed, to the best of my knowledge the three dimensional Thirring model 
may be the simplest fermionic model requiring a numerical solution, and thus de¬ 
serves study by lattice theorists interested in simulating dynamical lattice fermions. 

Lattice studies of the Thirring model have been made by two groups; I will 
briefly summarise xrur resultsS details of the other group’s approach are given by 
the next speaker. Ej The lattice formulations we have used differ in detail. I will 
suppress technical details, though it is still not clear how important they might turn 
out to be; one might hope that the number of distinct universality classes for fixed 
Nf consistent with the breaking of the symmetry (|2^) is small. We have studied 
lattice models corresponding to Nf = 2, i and 6 on volumes 8 ^, 12^ and 16^, with 
bare fermion masses ma = 0.05,..., 0.01 (a ~ is the lattice spacing). 

First let’s discuss the evidence for spontaneous chiral symmetry breaking. We 
have measured E = (ipip) various m and g'^, and fitted our data to the equation 
of state (EOS) 


The exponent 6 is identical to that of (M) - the form (|29| ) makes the implkit 
assumption that 6 — 1/P = 1 , which is consistent with the ladder approximationcS. 
To explore it is convenient to display the data in the form of a Fisher plot, 
ie., versus m/E. For the mean field exponents (W), this plot yields curves of 
constant as evenly spaced parallel straight lines, intercepting the vertical axis 
for g"^ > gl, the horizontal axis for g"^ < g^, and passing through the origin at the 
critical coupling. Any deviation signals a departure from the predictions of mean 
field theory. In Figs. 7, 8 , and 9 we show Fisher plots respectively for Nf = 2,4 and 
6 . The main features to note are that the Nf = 2,4 data show different curvature, 
and that the Nf = 6 data appear to be accumulating in the sub-critical region. The 
results of a fit to (p^), which include a finite volume scaling analysis for Nf = 2, 
are summarised in Tab. |^. 

We conclude that the Nf = 2 and Nf = 4 models do have fixed points, and that 
even if we relax the assumptions implied in (p9|) they fall in distinct universality 
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Figure 11: 


classes; the data collected so far support 4 < Nfc < 6, 

Next consider susceptibilities, which are equivalent to spatially integrated two- 
point functions. By analogy with ferromagnetic models we dehne a “longitudinal” 
(ie. scalar) susceptibility xi and a “transverse” (ie. pseudoscalar) susceptibility Xt 
as follows: 


XI = ; Xt = ^('075'0(O)V'75V'(a;))- (30) 

X X 


The longitudinal susceptibility has contributions from diagrams with both con¬ 
nected (□) and disconnected (o) fermion lines (ie. flavor non-singlet and singlet 
channels respectively), plotted in Fig. 10. The relative importance of the discon¬ 
nected piece, which is entirely due to sea quark loops, rises significantly in the 
broken symmetry phase. Fig. 11 shows that in the symmetric phase xi (°) and Xt 
(□) are approximately equal, as required by chiral symmetry, whereas in the broken 
phase Xt ^ Xh since Xt cx and we expect the pion to be a Goldstone boson 

in the chiral limit. The ratio xitxt = is another interesting quantity to 

examine. In general, for fixed g^, R varies strongly as a function of to as m —> 0, 
tending in the chiral limit to zero in the broken phase and to a constant ~ 0(1) in 
the symmetric phase. Exactly at the critical-coupling, however, it follows from the 


EOS (29) and the chiral Ward identity that 


= 


1 


(31) 


independently of to. Thus numerical measurements of R give an independent esti¬ 
mate of the exponent <5. In Fig. 12 we show data for Nf = 2 in the vicinity of the 
critical coupling, for 1/g^ = 1.9 (o), 1.92 (□), and 2.0 (o), together with the value of 
I from Tab. |^. The data suggest that R may becmie independent of to at criticality 
once finite volume effects are taken into accounts. Whilst this result is preliminary, 
it is interesting because measuring R may be the most promising method of testing 
the equivalence of the Thirring model with QEDg in future simulations. 

Finally I present results on the mass spectrum for Nf = 2, both for the physical 
fermion / and for scalar and pseudoscalar // bound states. A 16^ lattice is really too 
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Figure 13: 


small in the timelike direction for any great precision; however Fig. 13 demonstrates 
that the level ordering of the states varies sharply across the transition. There is 
no evidence for the absence of light composite states iri—the symmetric phase, as 
predicted for the case of a conformal phase transition E3. One further result is 
that in the vector channels in the symmetric phase there are light states in both 
“vector” i'ip'yfj.ip) and “axial pziector” ('i/' 7 /i 75 V') channels: the latter is not present in 
the large- approximationE3. 

So, the Thirring model is a good theoretical laboratory for non-perturbative 
dynamical symmetry breaking. In future studies it will be interesting to refine 
the spectroscopy, perhaps to explore the effects of a strongly-interacting continuum 
above the // threshold. More effort is still needed to pin down the value of Nfc, 
and to see if the Thirring model and QEDg are equivalent descriptions of the same 
critical dynamics. 

4 The Nambu — Jona-Lasinio Model for d = 4 

After all these nice results, it is natural to demand what happens in the physical 
dimension d = 4. Once again we revert to the GN model in the large-A/ expansion. 
As d —> 4, the expansion loses the property of renormalisability, essentially because 
the “vacuum polarisation” contribution to the auxiliary two-point function gets a 
subleading In A divergence. This must be compensated by the introduction of a 
kinetic term (cl^cr)^ in the renormalised Lagrangian. If instead we choose to work 
with a finite cutoff A, then it is possible to make do without introducing the scalar 
kinetic term, at the cost of a residual logarithmic dependence on the UV scale. A 
moderate separation of UV and IR scales is possible - the model is thus an effective 
theory, making continuum-like predictions but with 0(ln A) scaling violations. The 
UV scale can be interpreted as that at which new physical information, in the 
form of a more fundamental theory, must be supplied. The technical name for this 
scenario is triviality | - in the continuum limit A —> oo all interaction strengths 

*^Actually the model remains trivial even once promoted to a Higgs-Yukawa theory with the 
addition of (O^cr)^ and operators. 
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vanish, as we shall see below. In RG language, for d = 4 it is now the Gaussian 
fixed point which dominates the low-energy physics, and different models differ only 
by logarithmic corrections from free field theory. 

The model we will consider in detail is the Nambu - Jona-Lasinio (NJL) modelQ: 


Cnjl = - 


2Nf 


p'4’ip) {,'4’ip'lb'Tpq'^iq) 


(32) 


where p, q run over an internal isospace. In the chiral limit Cmjl has an SU(2)^ ® 
SU(2)^ axial symmetry, which at strong coupling is dynamically broken to SU(2)y 
by the generation of a condensate {'ip'4’) ■ This is similar to the pattern of symmetry 
breaking in the Higgs sector of the Stanciard Model, and has led to suggestions 
that the physical Higgs is a // composite Q. The field theory conventionally used 
to describe the Higgs sector is, of course, the 0(4) linear sigma model, describing 
elementary scalars with a four-point interaction. Both theories are trivial, and once 
the extra operators corresponding to scaling violations are takenjnto account, then 
each has equivalent predictive power for the Standard Models. That said, it is 
interesting to ask whether there is a distinction between theories such as the sigma 
model in which triviality is manifested through elementary scalars, and theories like 
the NJL model in which fermions are elementary but scalars are composite. 

One way of understanding the distinction is via the wavefunction renormalisa¬ 
tion constant Z, defined as the coefficient of I/fc^ in the scalar propagator in the 
IR limit 0. For the sigma model, Z is perturbatively close to one, whereas 

for the NJL model in the large-iVy expansion, for d < 4 Z oc (eg. (||)), and 

for d = 4Zocl/|lnI]|. In either case Z vanishes in the continuum limit S —> 0, 
with S measured in cutoff units: this is the compositeness condition llZl. Another 
difference lies in the relative ordering of scales. The interaction strength in the 
scalar sector is expressed in terms of the pion decay constant /tt, which must be 
finite for interactions to occur. For the sigma model, 

^^0(1), (33) 


however the ratio between the order parameter S and the physical Higgs mass 
is 


E I 
Ma 


lnE|i, 


(34) 


where gn is the renormalised (j)^ coupling constant, which vanishes logarithmically 
in the continuuniJimit. For the NJL model, on the other hand, to leading order in 
l/Nf we have 00 


S 


lnE|^ 



(35) 


All these predictions follow from the PGAG relation Z = E^. In both cases, 
triviality is realised by the divergence of /,r in units of the physical scale in the 
continuum limit: 


A 

M, 


lnE|i 


(36) 
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Table 4: The correction exponents q for d = A 



sigma 

NJL 

‘Thirring’ 0 

‘GN’B 

9i 

0 . 4 , 0.5 

0 

0.36(11) 

0.016(11) 

92 

1 

-1 

0.485(17) 

-0.553(18) 


The different triviality scenarios also predict different forms for the equation of 
state. Recall that in d = 4 the critical exponents take mean field values (|^). Due 
to triviality the terms in the EOS (^ ) are modified by logarithmic corrections, to 
give the form 

to = a("4-4 Vi ■ (37) 


IlnEKi 


|lnE|92 


In the sigma model, the powers qi and <72 can be found analytically by using 
RG to evolve the coupling from a generic value into the perturbative regime 
this works because even the perturbative EOS predicts symmetry breaking. For 
the NJL model, by contrast, symmetry breaking only ocxurs for large coupling 
strengths, where the only guide is the large-iVy- approachHil. The predictions are 
shown in Tab. (the two RG values quoted for qi correspond to U(l) and SU(2) 
symmetries respectively). 

l_Car^ctions to the mean field EOS have been analysed in several recent stud- 
iesE3^^, which use two different lattice implementations of the four-fermi inter¬ 
action. The first, which in d = 3 we have used to define the_lattice GN model, 
formulates the auxiliary helds a on the sites of the dual latticed, so that each •ifi'ip 
pair interacts with all other sites sharing a hypercube. This hasJaeen used in a 
study of the Nf = 2 NJL model with SU(2)^ O SU(2)^ symmetryci. The second, 
which in d =^2. defines the lattice Thirring model, formulates the auxiliaries on 
lattice links cZl; the Nf = 4 model with U(l) chiral symmetry was then studied^. 
The fitted values of the q exponents are also given in Tab. We should com¬ 
ment that logarithmic corrections are hard to isolate numerically, and that results 
are sensitive to assumptions about the size of the scaling region about the fixed 
point O. Nonetheless, the values found do suggest that the ‘GN’ and ‘Thirring’ 
formulations have qualitatively different scaling violations, the former resembling 
the large-predictions and the latter the sigma model. How stable the quoted g’s 
are, and whether the lattice-regularised models fall into distinct universality-elasses 
(the notion is still valid, since corrections to scaling can also be universal^), are 
issues for further study. 

Further support for the distinct nature of models with composite scalars comes 
from a study of finite volume corrections when the ‘GN’ simulation is run directly 
in the chiral limit m = OcHI. On a finite system in this limit the order parameter E 
is strictly zero; we must instead monitor the quantity |4>| defined by 


|$| = ( yitr$t$) ; $ = 4X!('^(^) 


(38) 


2 / ’ V- 

X 

where a, if are the auxiliary fields. In Fig. 14 we plot A = |$| — Eg versus 
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Figure 14: Finite volume correction for the NJL model. The vertical line shows the estimate for 

Vgl 


where S 
limit. Fc 
the form 


is the order parameter from the EOS extrapolated to the chiral 
.a theory with fundamental scalars, the correction is predicted to be of 


A cx 


1 


(39) 


where L is the linear size of the system {L = 16 in this case). For a composite 
scalar, however, in the large-Vj approximation the prediction becomes 


A cx 


1 

r^SollnSol' 


(40) 


As shown in the figure, (|40|) provides a much more convincing fit to the data. 


5 Conclusions 

To conclude, let me briefly summarise what each of the models we have looked at 
has taught us: 

• The Gross-Neveu models with discrete and continuous chiral symmetries in 
d < A appear under good analytical and numerical control. They provide 
a paradigm for a non-perturbative fixed-point theory containing elementary 
fermion and composite scalar degrees of freedom, and a non-vanishing inter¬ 
action at short distances. Two further applications of these models, which I 
have not had time to discuss, is that firstly unlike gauge theories they per¬ 
mit Monte Carlo simulation at non-vanishing chemical potential, and hence 


16 






non-vanishing baryon number density. This has been applied in d = 3 as-a. 
model effective theory of hadronic physicscj. In d = 4 simulations of QCDH 
directly in the chiral limit are possible by including a gauge-invariant four- 
fermi interaction in the bare Lagrangian; the resulting model lies in the sapie 
universality class and hence contains the same IR physics as the true theoryEil. 

• The Thirring model in d < 4 has a fixed point for small Nf not described by 
the 1/Nf expansion. Genuine non-perturbative techniques need to be applied, 
and the model deserves to be used as a testbed for numerical lattice studies. 
The issue of its equivalence with QEDg is intriguing and needs to be resolved. 

• There is lots more to learn about the nature of triviality in d = 4. So far in 
looking at corrections to the equation of state we have pursued rather a “con¬ 
densed matter” approach, and argued that composite models have a distinct 
form of scaling violation as compared to the standard ferromagnetic picture. 
It would be interesting to try to extract particle physics from these models, for 
instance, by asking whether they predict new triviality bounds on the Stan¬ 
dard Model Higgs mass. Finally, there remains the fascinating possibility of 
actually finding new non-trivial fixed points in d = 4 by expanding the space of 
couplpgs, say by introducing a gauge interaction such as in strongly-coupled 
QEDEi 
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